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We consider the model of superconducting pairing with the energy gap function 
which is odd over k—kp. In this case superconductivity is possible even in the presence 
of an arbitrarily large point-like repulsion between electrons, which is attractive for 
the theory of high-temperature superconductors. We discuss mainly a model pairing 
interaction for which the BCS equations can be solved exactly, allowing the complete 
analysis of the interplay of the usual ("even") and "odd" pairing, depending on the 
coupling constant. We also show that the normal impurities (disorder) lead to a rapid 
degradation of the "odd" pairing, which is even more rapid than in case of magnetic 
impurities in traditional superconductors. 
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I. INTRODUCTION 



Traditional BCS theory of superconductivity is based upon an assumption of 
the existence of an attractive interaction of electrons with opposite spins close to 
the Fermi level. Usually it is assumed that this attractive interaction overcomes the 
Coulomb repulsion of electrons at least in some part of the phase space, which is con- 
sidered to be the necessary condition for superconducting state at low temperatures. 
Naturally, the strong repulsion of electrons characteristic to Hubbard-like models, 
which are extensively used to describe the electronic properties of metallic oxides, is 
a factor opposing superconductivity. From this point of view it is interesting to con- 
sider models, in which the influence of such repulsion is strongly suppressed or even 
absent. A model of this kind was proposed in a recent paper by Mila and Abrahams 
0. It is based on the assumption that the BCS gap A(k,u) depends only on \k\, 
or more precisely on the quasi-particle energy £ fc = — kp), measured from the 

Fermi level and is an odd function of this variable. In this case superconductivity 
becomes possible even for the arbitrarily strong point-like repulsion between elec- 
trons. Physically it is obvious that such a state may realize in case of strong enough 
repulsion, when the usual ("even") superconductivity is suppressed. References to 
the earlier works on the "odd" -pairing can be found in Ref. [[3]. 

The aim of our paper is more detailed than in Ref. || study of the "odd" pairing 
and its relation to the usual "even" case, using the simplest weak coupling ap- 
proximation and a model pairing interaction which allows the exact solution of BCS 
equations. This model allows the detailed analytic analysis and comparison with the 
results of Ref. || which were obtained numerically for more "realistic" interaction. 
Besides that we consider the influence of normal impurities on the "odd" pairing. 
This influence is unusually strong || and superconductivity is completely suppressed 
even faster than in the case of magnetic impurities in traditional superconductors. 
We solve this problem both using the model pairing interaction and numerically in 
the case of "realistic" interaction proposed in Ref. Despite the obvious attrac- 
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tiveness of the odd pairing for the analysis of high-temperature superconductivity 
in oxides its strong instability towards disordering makes it unlikely candidate for 
pairing in cuprates. 



II. EQUATIONS FOR THE GAP AND 
TRANSITION TEMPERATURE 

The model is based on the demonstration of the fact || that the weak coupling 
BCS gap equation ffl: 



A(0 = -N(0) 7<%V&0 , A(e ' } th [ >/t a+ fM ) (i) 



can acquire the nontrivial solution A(£) = — A(— £) (i.e. odd over k — kp, £ = 
vp{k — hp)) in case of the presence in V (£,£') of a strong enough attractive part 
even if very strong (infinite) point-like repulsion is also present. It is easily seen |2| 
that for the odd A(£) the repulsive part of interaction in ([!]) just drops out while the 
attractive part ^(£,£0 may lead to the pairing with some unusual properties: the 
gap A(£) is zero at the Fermi surface which leads to the gapless superconductivity. 
We must stress that in this case we have isotropic pairing and the gap is zero 
everywhere at the Fermi surface which is different from the case of anisotropic, e.g. 
d-wave pairing 0. 

Thus, in the following we assume that the interaction kernel in Eq.(|l]) consists 
of two parts: V{£, £') = K(£, £') + V 2 {£, £'), where 

[ U > for |£|, |£'| < E F 
V^X) = \ (2) 
[0 for |£|,|£'|>^ 

— is the point-like repulsion of electrons, while ^(£, £') — is an effective pairing in- 
teraction (attraction), which is nonzero for |£|, |£'| < uo c and |£ — £'| < uo c (the last 
condition is of the major importance) where u c Ep — plays the role of char- 
acteristic frequency of Bosons, responsible for the pairing. The pairing "potential" 
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V2(£, O can be modelled by various functions, e.g. by a step-function as in Ref. [0. 
In Ref. the major attention was to the following form of this interaction: 

\ -v{\i-i'\/uo c y 2 i z for leue'ue-ei <E F 

V 2 (Z,0 = \ (3) 

[o for ieueiorie-ei>^ 

which was justified by the wish to get the following tunneling density of states of 
a superconductor: Nt(E) ~ \E\ 2 for E — > 0. The integral gap equation (|l|) with 
such V^(f ,£') was solved in Ref. |2| numerically, and some other results qualitatively 
similar to those observed in high-temperature superconducting oxides were obtained. 
Here we shall concentrate mainly on the following model interaction ||: 

\ -VcosU^-) for |£U£'U£-£'|<w c 
[0 for |£|,|£'|or|£-£'|>u; c 

The major advantage of such a choice is that now we can reduce the integral gap 
equation to a transcendental one and solve it easily. The choice of Eq. (|]) is not 
in this sense unique, quite a number of similar "potentials" may be proposed. We 
can use for example ^(£,£0 ~ c/i(£ — £') or (£ — f') 2 . However, in case of Eq.(|J) 
we obtain the results which in some sense are closest to those obtained with the 
"realistic" choice of (|3]) (e.g. for the density of states). The majority of qualitative 
conclusions below do not depend at all on the form of the model "potential" . The 
importance of the model interaction of Eq.(|4]) is obviously related to the fact that 
any "potential" which is an even function of £ — £' on the interval from — u c to u c , 
can be represented by the Fourier expansion over cosines. In this sense our analysis 
forms the basis for rather general study. 

We should like to note, that the real choice of the pairing interaction V(^, £') can 
be made with the use of dielectric function formalism of superconductivity theory 
I^U, where rather general expressions for the interaction kernel of the BCS theory 
are obtained. Unfortunately, it is unclear how we can get any nontrivial dependence 
of this kernel upon |£ — because in this formalism we typically get the dependence 
ofV(£,?) on |f | and separately |,|. 



Superconducting transition temperature is determined by the linearized equa- 
tion: 

A(0 = -N(0) Jd?V(te^^p-th (jfr) (5) 

— oo ^ ' 

Using (0) and (||) we get: 



for |f | < u c , while for oj c < |f | < E F 



(7) 



where we have defined the usual dimensionless constants of pairing and repulsive 
interactions: g = N(0)V, fi = N(0)U. 

The general solution of Eqs.(H), (H) takes the form: 



A(0 



A C co S (f£)+A sS m(f£)+A for |£| < 



(8) 



A for lu c < |f | < Ep 

where A c , A s , A are determined by the following algebraic system of equations: 

A c = gF c A c + gFA 
A = -fiF A c - y. W'A 



(9) 



A s = gF s A s 



(10) 



where: 



F s 
F 
W 

w 



_L 

2T C 



f dices* (f£) \th{ 





U! c 







2T C 



fd£±th{£r) 
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We can see that the equation (0), determining T c for the "odd" pairing is de- 
coupled from the equations (0), determining T c for the "even" case. The repulsive 
part of interaction influences only on the "even" pairing and from Eq. (0) we get the 
following equation for T c : 

F 2 

l=gFc- 9^ (12) 

which may be rewritten as: 

l=gF c - n*W + /2*g(F c W - F 2 ) (13) 
where we have introduced the Coulomb pseudopotential: 



l + fi(W'-W) 

where in the weak-coupling region, when T c -C w c the value of W — W reduces to 
the usual In— . 

Transition temperature for the "odd" state is determined by: 



1=^ = 9/^1^1-^1^1 (14) 







In the Appendix we derive Eqs.(|l2D and ([14]) from the usual picture of Cooper 
instability of the normal state. 



In Fig.l we show the results of numerical solution of Eqs.([L2j) and ([14]) for dif- 
ferent values of coupling constants g and //. For weak repulsive interaction "even" 
pairing dominates and the appropriate transition temperature is higher than the 
transition temperature for "odd" pairing. As repulsion grows and for larger values 
of g the "odd" pairing overcomes the usual one. Note that for the model interaction 
of Eq. (f|) we have the critical value of the pairing constant and the "odd" pairing 
appears only for g > g c w 1.213. Thus we actually move from the region of weak 
coupling for which the BCS equations are derived. In this sense , the results shown 
in Fig.l for the region of large coupling constants are more or less of illustrative 
nature only. In particular, this concerns the practically linear growth of T c with 



the growth of g, which is obvious from Fig.l and is apparently due to uncontrolled 
extrapolation of BCS equations to the region of strong coupling and large g. In fact 
here we must use much more elaborate analysis in the spirit of Ref . , where the 
smooth crossover from "large" Cooper pairs of the weak coupling theory to "com- 
pact" Bosons of the strong coupling region was first demonstrated. It is known, 
that in this case the growth of g leads to the saturation of T c , which for the large 
g region is determined by the well known expression for the temperature of Bose 
condensation in an ideal gas of Bosons with practically no dependence on g. 

The critical value of the coupling constant g c for the "odd" pairing is formally 
absent if we use the pairing potential of Eq.(|3|), which is obviously due to its di- 
vergence for |£ — £'| — > 0. The appropriate dependence of T c on g = N(0)V can be 
obtained by direct numerical solution of Eq.fl5|) with potential given in Eq.(^) and 
is shown on the inset in Fig.l. At the same time it is clear that in this case also the 
"odd" pairing starts to dominate over "even" only in case of strong enough repul- 
sion. Note, that the analysis of strong repulsion within the framework of BCS model 
is again rather doubtful, because it takes into account only the simplest Fock correc- 
tion due to electron-electron interaction. It is clear then that we cannot hope to get 
consider the limit of \x — > oo rigorously within this approach. However, the above 
analysis of BCS equation apparently produces the correct qualitative description of 
the crossover from traditional "even" to the "odd" pairing. 

Consider now the temperature dependence of the gap function for the case of 
the "odd" pairing for the model interaction of Eq.(^). According to Eq.(§) for the 
"odd" case we have: 



and the temperature dependence of A (T) is determined by the following equation 
which is easily obtained from Eq. ([!]): 




(15) 







for |£| > u c 
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,x th 1 1 



2T 



£' 2 + A§(T)sin 



2 (tlSL 

2 w c 



(16) 



Solutions of this equation for the number of values of the pairing constant g are 
shown in Fig. 2. The temperature dependence of A (T) resemble that of the BCS 
theory but is not identical to it. In particular, for the large coupling constants 
g 3> 9c the value of 2Aq(T = 0)/T c ~ 5 with some tendency to become smaller as g 
diminishes. 

The tunneling density of states is easily calculated in a usual way ||. Using (|15|), 
we get: 

ME) I r^ggf ■ , ,~\ for l e ' <LJc 

for to r < lei 



iVo(O) 



where e is determined from the equation e 2 + Ag(T)sm 2 (f ^7) = -S 2 - The appro- 
priate dependences for different temperatures are shown in Fig. 3. Density of states 
is always gapless, "pseudogap" is smeared as temperature grows, at the same time 
the positions of peaks are only weakly dependent on temperature. These results are 
qualitatively close to those obtained in Ref. for the case of interaction given by 
Eq. (|3]) and can be compared with the known anomalies of the gap-behavior in high- 
temperature superconductors. If we define the ratio ^ using the peaks positions of 

J- c 

the tunneling density of states we find ~ 6. 

J- c 

III. THE INFLUENCE OF NORMAL IMPURITIES 

The major interest is to study the influence of normal (nonmagnetic) impurities 
on the "odd" pairing. It is well known |l]|| that such disorder practically does 
not influence traditional "even" pairing. In our case the equations for normal and 
anomalous Green's functions have the usual form ||, which is valid in case of weak 
impurity scattering: 



S 



2 



F(u£) = =^4 (18) 



where to = (2n + l)7rT, 



do 



T f AC u 

oo = to / at,— 
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A(0 = A(0 + - 1 d£ = A(0 (19) 

Here 7 = 7rcV 2 iV(0) — is the scattering rate of electrons due to point-like impurities 
with potential Vq, chaotically distributed with spatial concentration c. The integral 
in the second equation in fll9|) is actually zero due to the odd nature of A(£), so 
that gap function renormalization due to impurity scattering is absent. This is the 
main reason of strong impurity effects in case of the "odd" -pairing. Note that the 
similar situation is characteristic for anisotropic, e.g. <i-wave pairing P,|10|. 
Gap equation takes now the following form: 

A(0 = ~iV(0)T£ Jd^V^ f^ 2( (20) 

Close to T c this equation can be linearized: 

a(o = -iv(o)r£ / dev^a^rU (21) 

Here Co = 00 + 'ysignuo. 

The sum over Matsubara frequencies in (|2T|) can be calculated in the usual way 
transforming into integral in the complex frequency plane. As a result, the linearized 

gap equation can be written in several equivalent forms: 

00 00 

A(O = -iV(0) J deV 2 (£,0 / ^th(^ReG R (-uje)ImG R (ue)A(0 (22) 

—00 —00 

where G R (u£) = {00 — £ + — is the retarded Green's function of a normal metal 

with impurities. In another way we may write: 

00 00 . . 

A( = -N(0) / ifVitt, O / ^tfc ^ _X_ A(f ) (23 ) 



Equation similar to (p2|) was recently obtained in Ref.( fLOp for superconductors 
with anisotropic <i-wave pairing, in the following we use it in the form of Eq. fl2"3|) . 

For the model interaction of Eq.([|) the gap again takes the form given in Eq.(|i~5|) 
and T c -equation directly follows from Eq. (|23|) : 



1 =»/t 4ta U^yT t H-iirJ^ (24) 

In Fig. 4 we show the dependences of T c on 7 for a number of values of the pairing 
coupling constant g which were obtained solving Eq.(^). It is seen that the nor- 
mal impurity scattering strongly suppresses the "odd" pairing. Superconductivity is 
completely destroyed for 7 ~ T c o, where T c o — is the transition temperature in the 
absence of scattering (7 — > 0), defined from Eq. (|1"4"D . The destruction of supercon- 
ductivity occurs here even faster than in case of introduction of magnetic impurities 
into traditional superconductor |Tl[. In particular this is demonstrated by the fast 



disappearance of superconductivity region on the "phase diagram" of Fig.4 as g — > g c 
and also by the absence of universal dependence of T c (7), characteristic for the case 
of magnetic impurities. 

In case of the model interaction given in Eq.@ the T c dependence on 7 can 
be obtained by direct numerical solution of the linearized equation (fffi). To cal- 
culate the minimal characteristic value, determining the coupling constant g for a 
given T we have used the traces and Kellog's methods fll2"f . Integrals containing 
~ |£ — £'|~ 2 / 3 were calculated by methods which were proposed for singular integrals 
|j~3|l , allowing to estimate such integrals with an accuracy of the order of usual Gaus- 



sian quadratures. The procedure of calculation of minimal characteristic values was 
rather sensitive to the accuracy of calculation of symmetrized kernels. Satisfactory 
results were obtained using the representation of these kernels via hypergeometric 
functions, which were calculated by summing the appropriate generating series with 
a given accuracy. Kellog's method while being pretty fast in comparison to method 
of traces has a tendency to instability in the region of small coupling constants. 
The results obtained for T c (7) are shown in Fig. 5. It is obvious that the qualitative 

10 



picture obtained from the simplified model interaction above is conserved for this 
more "realistic" interaction. 

The critical scattering j c , corresponding to the total destruction of superconduc- 
tivity (T c (7 — > 7 C ) — >• 0), is determined by the following integral equation: 



A(0 = -N(0) J dgV^^-Larctg A(£') 



(25) 



which for the interaction given in Eq.(01) reduces to: 



iJt^KWfe) (26) 







For g ~ g c from here we can get the dependence like 7 C ~ (g — g c ) — > 0, which 
describes the disappearance of superconductivity region in Fig.4. For g ^ g c (strong 
coupling) we get the universal result: 7c/T c0 = 4/tt m 1.273. In fact this result as 
well as the dependence of T c (j) for g ^> g c are independent from the choice of model 
potential V2 (£,£')■ Concerning the universality of the ratio •fc/Tco, it is eas Y to see 
that it follows from the equivalent form of Eq.(|l4"|) for T c0 and Eq.(^) for 7 C in case 
of T c0 u) c and 7 C 3> u) c (i.e. for g ^> g c ): 



A(0 = -^ / rf£V 2 (£,OA(0 (27) 



-co 



a(o = -^ /^v 2 (e,r)A(r) (28) 

7T7r •/ 



-co 



The equivalence of these equations allows us to equate the coefficients of integrals 
which immediately gives the above result for the ratio of 7 c /T c o. Analogously it is 
easily seen that Eq. (p3|) for T c (j) ^> u c reduces to: 



00 



a(o = -^rKf ( 7?rrT ) I dev a (e,OA(o <- X) > 



4T C ( T )' V T e(7) 

where 

/(a 



du 1 x 



vr c/i 2 (o;/2)x 2 + ^ 
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Accordingly, comparing Eqs. (ref32) and fl27|) we see that in the limit of strong 
coupling the dependence of T c on 7 is determined by the following universal equation: 

T f } 7 , = 1 (30) 
T cof {jfa) 

We must stress however, that these "strong coupling" results are of rather formal 
nature and should be seriously modified in the spirit of Ref . @ . 

Our results for j c for the case of the model interaction of Eq.(|J), and those 
obtained by numerical solution of Eq.fl2~5|) with "realistic" interaction of Eq.(||) are 
shown at the insets in Fig.4 and Fig.5. 

IV. CONCLUSION 

In conclusion let us formulate the main results. We have presented the simple 
model of the pairing interaction which allows the complete analysis of integral equa- 
tions of the BCS theory both for the usual "even" (over k — kp) and for the exotic 
"odd" pairing. It is shown that the "odd" pairing becomes preferable for sufficiently 
strong electron repulsion and, in general, for sufficiently strong pairing interaction. 
The last problem ( "strong coupling" ) deserves further studies of the crossover from 
Cooper pairs to compact Bosons. The "odd" pairing leads to the picture of gapless 
superconductivity and to other deviations from the traditional BCS theory such 
as the unusual evolution of the pseudogap in the density of states, the large ra- 
tio of 2A /T C etc., which are attractive from the point of view of the theory of 
high-temperature superconductors. 

At the same time normal impurities (disorder) strongly suppress the "odd" pair- 
ing. This effect is even stronger, than in the case of magnetic impurities in usual su- 
perconductors. Even in the "strong coupling" limit superconductivity is completely 
destroyed for 7 ~ T c q, while for the smaller pairing constants superconductivity 
region on the phase diagram rapidly shrinks. 

It is rather well known that high-temperature superconductors are unstable with 



respect to structural disordering |]14| . However, if we exclude special cases like the 
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impurities of Zn, superconductivity is suppressed in metallic oxides rather close to 
the disorder-induced metal-insulator transition, which is apparently due to Ander- 
son localization |TJ]]. Anderson transition takes place for 7 ~ Ep S> T c o, so that 
the "odd" pairing is destroyed well before. Apparently this makes the "odd" pairing 
rather unlikely explanation of high-temperature superconductivity in oxides. How- 
ever, we can not exclude the possibility that some effect in these systems can be 
explained by the fast suppression of the "odd" component of superconducting or- 
der parameter by disordering, while the "even" part which is stable to disordering 
remains. This problem deserves further studies. 

This work was supported by the Scientific Council on High- Temperature Super- 
conductivity under the Project N° 90135 of the State Research Program on Super- 
conductivity. It was also partly supported by the Soros Foundation grant awarded 
by the American Physical Society. 
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APPENDIX 



It is useful to derive equations determining T c both for "even" and "odd" pairings 
from the picture of normal state instability, i.e. as the points of the relevant Cooper 
instabilities. Consider the two-particle Green's function in Cooper channel which 
is shown diagrammatically in Fig. 6. It is convenient to analyze instability of the 
normal state associating it with the divergence of this Green's function, summed 
over Matsubara frequencies: 



for q = Q = 0. Consider again electron-electron interaction V (£,£') consisting of 
two parts defined in Eqs.(||) and (f|). Due to isotropy of the system $ pp /(00) may be 
represented as $(£,£'), which is defined by the equation: 



(31) 



oo 




(32) 



— oo 



where 




(33) 



Taking into account Eqs.(|2|) and (§) we have: 




(34) 



f° r < u c and, accordingly: 




(35) 
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for |e| or |£'|>u; c , \Z\,\?\<E F . 
Let us define the functions: 



no 



(36) 



-U) C 

-E F 



for which, using Eqs.( 



we obtain the following system of equations: 



/(0 = Z(£)-nW'f(t)+gFf e (£) 

/ c (0 = Z(Oco S (f£)-//F/(0 + <?F c / c (0 

/.(O = z(e) S m(|i;)+^F s / s (e) 



(37) 



where we have used notations introduced in (|iT|). 

Now we can see that the "odd" and "even" equations has decoupled. The "odd" 
pairing is connected with the divergence of the function / s (£), and the corresponding 
instability condition has the form of 1 = gF s , which coincides with Eq. fll4|) . The first 
two equations in (|3T[) determine the instability towards "even" pairing. It is easy to 
see that 

and instability condition takes the following form: 

F 2 

1 = gF c - gn 



1 + fiW 



(38) 



which coincides with Eq.( 
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Figure Captions: 



Fig.l. Dependence of T c0 on the pairing coupling constant g = N(0)V in case of 
the "even" (lines) and "odd" (dashed lines) pairings. 

1 — n = 0; 2 — n = 1; 3 — ji = 10. In numerical procedures it was assumed that 
E F /u c = 50. 

At the inset — analogous dependence for "realistic" interaction of Eq.(3). 

Fig. 2. Temperature dependence of the gap A (T) in case of the "odd" pairing 
for the number of values of pairing coupling constant. 
l—g = 1.22; 2—1.5; 3—2.0; 4—3.0; 5—4.0; 6—5.0. 

Fig. 3. Density of states in the model of the "odd" pairing for a number of char- 
acteristic temperatures. 

1— T/T c = 0; 2—0.6; 3—0.8; 4—0.9; 5—0.99. 
Pairing coupling constant was assumed to be g = 3. 

Fig. 4. Dependence of T c for the "odd" pairing on the impurity scattering rate 7 
for different values of the coupling constant g: 

l—g = 1.22; 2—1.24; 3—1.30; 4—1.50; 5—2.0; 6—5.0; 7—10.0. 

At the inset — dependence of the critical scattering rate on the pairing coupling 
constant. 

Fig. 5. Dependence of T c for the "odd" pairing on the impurity scattering rate 7 
for different values of the pairing coupling constant g in the model with "realistic" 
interaction of Eq.(3): 

l—g = 0.17; 2—0.25; 3—0.72; 4—1.15; 5—6.41. 

At the inset — dependence of the critical scattering rate on the pairing coupling 
constant 

Fig. 6. Two-particle Green's function in Cooper channel. 
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